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Abstract 

We continue to develop further the bag-type Double Phase Transition Model (DPTM) 
for transformation of Quark-Gluon Plasma {QGP) to normal hadronic matter (H-phase). 
The model is based on the assumed existence of an intermxdmte Q phase composed by 
massive constituent quarks and pions (as Goldstone bosons) El’H. 

In the present paper the problem of entropy S and baryon number Mb conservation 
in phase transitions from deconfined phases {QGP and Q) to hadronic matter H is consid¬ 
ered. It is shown that standard construction of both first order phase transitions, H ^ Q 
as well as Q QGP implies a discontinuous structure of entropy per baryon S/Mb when 
crossing phase boundary; this results in impossibility of equilibrium transition from QGP 
to hadron gas. 

We follow the way suggested recently lid’ by H.Satz et al. for the same problem 
concerning direct transition H QGP. They proposed a modification of bag pressure 
parameter Bqgp by making it dependent on system temperature T and baryon chemical 
potential fi; this modification has been demonstrated to be sufficient to provide S/Mb 
conservation. 

Here we show that within DPTM such a modification turns out to be necessary 
and sufficient for bag pressure Bq in the Q phase only. The DPTM modified in such 
a way is shown to satisfy equilibrium Gibbs criteria for phase transitions. Location of 
phase boundaries in fi — T plane has been demonstrated to be changed but slightly; the 
modification tells mainly on baryon number density within Q phase. Two alternative 
descriptions of nucleon-nucleon interaction -the Hard Gore Model and the Mean Field 
Approximation - have been tested; the results for both cases appeared to be similar. All 
the results are shown to be stable against rather broad variations of model parameters. 
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1 Introduction 


According to the fundamental QCD predictions strongly interacting matter has to exist in 
(at least) two different phases El : 

- Quark-Gluon Plasma (QGP), i.e. gas of deconfined and massless quarks q, antiquarks q and 
gluons g, - at super high temperature T and energy density e, and 
-Hadron Gas , or iJ-phase at low T and/or e. 

QGD lattice calculations □ confirm this concept indicating to the first order phase transi¬ 
tion at T 200 MeV. Well above this temperature matter behaviour is easily described by QCD 
perturbation theory. Close to the critical conditions perturbation theory fails. However critical 
behaviour of strongly interacting matter is very important to analyze the phase transition condi¬ 
tions and the QGP signals. To study this problem one has to use some phenomenological models. 
Thermo^namic approach with two-phase matter Equations of State (EOS) has been widely used 
recently 0 and gave a lot of interesting results and experimental predictions. However there still 
exist some ambiguous points connected with phase transition description. One of them is discussed 
here. _ 


We consider, following the common wayQ, bag type EOS for QGP and-usual nonrelativistic 
EOS for H phase (taking into account hadron interaction). Such type modelsB are known to result 
in first order phase transition H <-!■ QGP at some critical temperature T^, if none intermediate 
phase is taken into account (to be called later on Single Phase Transition Model, or SPTM). 
However, in distinction to the great majority of researchers considering only direct phase transition 
H <-> QGP, and thus assuming that, both, quark-deconfinement and chiral symmetry restoration 
have to occur simultaneously, we have worked outEl’O the Double Phase Transition Model (DPTM) 


based on the assumption of an intermediate phase Q formed by 
pions (as necessary Goldstone bosons). In the form proposed in 
that such a Q phase could exist. 

The idea of possibility of two phase transitions goes back to E. ShuryaklD 


nfined constituent quarks and 
it was actually demonstrated 


who put forward 


various arguments to point out that temperatures of quark deconfinement Td and chiral symmetry 


restoration Tch may not coiui 
by other works within QCDIS 


ie, and Td should be less than Tch- Later this idea was supported 
i. It leads to the conclusion that there may exist some temperature 
interval Td < T < Tch where quarks are liberated off individual hadrons but still possess non-zero 
mass. Such objects are well known from Additive Quark Model (AQM) and called constituent 
quarks, or valons. Xhey are necessary entities for satisfactory description of moderate energy 
hadron phenomena lj. Thus an intermediate phase of deconfined massive constituent quarks, 
Q-phase, may exist. 

Eirst attempt to investigate tl^ problem within the thermodynamic models with ba^ype 
EOS belongs to the Bielefeld group til’113; then this problem has been investigated in Refsli3“t3. 
A possibility for existence of the intermediate phase Q formed by deconfined constituent quarks 
and pions was in fact demonstrated. However the choice of the key model parameters based 
on the lattice calculation data for baryonless matter 0 resulted in the negative conclusion that 
H <-> QGP transition should proceed almost always directly, without any intermediate state, since 
in temperature T - chemical potential /i plane Q phase occupies only tiny petals and can hardly 
play any essential role in cecity. . j—. 

In our earlier worksQ’0 the problem has been reconsidered following the same ideologyllH” 
but with different physical approach to the choice of the bag parameters, since lattice approach (not 
securing pion description) does not seem to be a safe basis in the case of <-> Q transformation 
where pions play the decisive role. 

This resulted in conclusion that Q phase seemingly exists almost always (at least, for fi < 1). 
In p—T plane it occupies a corridor between H and QGP phases having the width AT = Tch — Td ~ 
50 MeV. The value of Td was found to be equal to some 150 MeV and its physical meaning (the 
highest temperature allowing for existence of hadrons, above T^they have to decay into constituent 
quarks) enables to identify it with the Hagedorn temperaturelij. Both phase transitions are of the 
first order. 

The qualitative stability of these results for varying values of parameters (within reasonable 
limits) has been demonstrated. More details are given in the section 2. 

The present paper is devoted to the problem of entropy and baryon number conservation 
when crossing phase transition boundaries. This problem appeared already in the common model 
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with single (first order) phase transition, SPTM. It has been shown til’ within SPTM that as a 
result of the transition H ^ QGP the entropy per baryon S/Nb is discontinuous across the phase 
transition boundary: 

{S/Nb)qgp > {S/Nb)h- 

This means that transition H <-> QGP is irreversible: it could not satisfy thermal and chemical 
equilibrium conditions, and, at the same time, fulfill baryon number and entropy conservation at 
the phase boundary. In particular, adiabatic transition QGP ^ H is impossible as it should be 
accompanied by the entropy decrease (since conservation of the baryon number is secured). 

Thus one has to choose between two possibilities: 

-either the transition from QGP into H-phase can not proceed as an adiabatic one under any 
conditions, 

- or EOS used are not fully correct. . 

The first possibility has been discussed recently Ell in connection with large value of specific 
entropy {S/Nb)h detected in experiments on high energy heavy ion collisions. It was used as an 
argument for assumption that QGP has been actually observed. 

The later possibility seems rather realistic due to uncertainties presenting in phenomeno¬ 
logical EOS. In particular, bag-type EOS used for QGP includes the key model parameter Bqqp 
(representing nonperturbative interactions of quarks and gluons with physical vacuum). In various 
works on SPTM its value was chosen rather arbitrary: Bqqp = 0.2 -h 0.5 GeV/fm^. Usually it is 
treated as constant, but there are no reasons for B not to depend on T and/or /x. 

Accordingly in E3 there was suggested a certain modification of Bqqp making it p and 
T dependent, B{p,T), instead of commonly used Bqgp = const. This enabled to restore the 
continuity of S/Mb and thus to solve the problem immediately. This modification of EOS was 
shown to change phase diagram of the system not considerably, while the transition QGP <-> H 
becomes reversible. 

In this paper we study the problem of specific entropy conservation within DPTM. There 
appears the same discontinuity of S/Mb when crossing both transition boundaries, thus both phase 
transitions appear to be irreversible. It is shown that the reversibility of both phase transitions 
can be restored by modifying the Q phase EOS alone, without changing Bqgp. 

The paper is organized as follows. In section 2 we remind basic features of the DPTM. Spe¬ 
cific entropy discontinuity and the method of its correction within SPTM is discussed in section 
3; in section 4 the same problem is discussed within DPTM; results of numerical calculations are 
presented in subsection. Section 5 represents summary and discussion. Some details of hydrody- 
namical description of equilibrium system evolution are given in Appendix. 


2 Basic features and main results of DPTM 

We use, following common way0“Ei, bag-type model EOS for QGP and Q phases: {p,gi, 
mi are pressure, degeneracy factors, masses of i-th type particles respectively for each j-th phase; 
j means hadronic i?,valonic Q, and QGP phases respectively). 
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Massless gluons {g) and quarks (q) participate in QGP; Q phase contains constituent quarks 
(m, d, s with ~ mj, —320 MeV and TOs — 512 MeV), and pions; p^Q being chemical potential of 
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constituent quarks, equal to that of corresponding current quarks [^iq = Hu = and fjLg is taken 
to be zero). 

The terms Bj in EOS of Q and QGP reflect effective interactions with vacuum. Bq, Bqqp 
are free bag parameters chosen according to their physical meaning: Bqqp is QCD vacuum energy, 
density knownE2l to be « 0.5 ^ 1 GeV/fm^, and Bq is estimated from low energy phenomenaral 
as 50 -h 100 MeV/fm^ (close t| 0 ,the B value of MIT-bag model). Note that this choice of B values 
(instead of that chosen inEdl“EJ) is crucial for appearance of the intermediate Q phase. 

For the hadronic phase H : 


PH{T,fi,V) = 




k^dk 

a/p ^2 
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The first term in (3) represents pion contribution; the summation (in the second term) is over all 
i — th type stable hadrons dominating in the H phase (tt, N, A and K were taken into account); 
rrii , Pi and pi are the corresponding masses, chemical potentials and degeneracy factors. The last 
term st^ds for account of nucleon-nucleon interactions in the form of Mean Field ApproximMion. 
(MFA) EiJ; we have used also Hard Core Model (HCM) for nucleon interaction description EJ’ El 
(where the form of interaction account in the H phase is more complicated). 

In thermodynamic equilibrium, according to Gibbs principle, the actually realized phase is 
that with the largest pressure at given p and T. Then at given value of nucleon chemical potential, 
p there are 3 possible transition temperatures: 

deconfinement transition (H <-*■ Q) at T^: pH(Td,p) = pplTd, p/3) 

direct transition {H ^ QGP) at T).: ph{Tc,p) = pqgp{Tc, p/3) 

chiral transition {Q ^ QGP) at Tch'- PqiTch, p/3) = pgapi/Tch, p/3) 

The case of coincidence of all three transitions corresponds to the triple point at 
Ph{T*, p) = pq{T*, p/3) = pqgp{T*, p/3) 

Q phase actually exists if for rising temperature the deconfinement of valons occurs first, 
prior to direct formation of the QGP phase, i.e. Td < T^; and, in the opposite direction (for 
decreasing temperature), if chiral transition from QGP to Q phase occurs prior to formation of H 
phase, i.e. Tch > Tc- Thus general condition for Q phase existence is: 

Td{p) < Tc{p) < Tch{p)- (4) 

In this case the direct transition does not occur. Otherwise DPTM reduces to a model with single 
phase transition and its results coincide with that of SPTM. This very case have been met in the 
papers llii“ El due to specific choice of the B parameters. 

The choice of ftdel parameters based on its physical meaning resulted in the quite opposite 
conclusions, namely B S: 

i) H ^ QGP transition proceeds almost exclusively via the Q phase, H ^ Q ^ QGP. 

ii) Deconfinement of valons PI ^ Q should occur at rather low energy density of nuclear matter 
« 0.3-;-0.4 GeV/fm^ r(only three times larger than energy density in a normal nucleus, as it was 
roughly estimated ini). 

iii) Temperature interval for Q phase, AT = Tch — Td, amounts to ^ 50 MeV (see Fig. 1). For 
baryonless matter, typically, Td « 140 MeV and Tch ~ 200 MeV. Thus Td coincides with the well 
known Hagedorn temperature (as it should be since at Td hadrons cease to exist and decay into 
constituent quarks), whhdi had been treated earlier as some approximation to direct transition 
H <-> QGP temperaturellj and now is shown to have actually independent physical meaning. 

These results turned out to be qualitatively stable against extended variations of model 
parameters and nucleon interaction description. 
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3 Restoration of specific entropy continuity in 
SPTM. 


As it has been said above, within SPTM the specihc entropy value turns out to be discon¬ 
tinuous when crossing the direct transition boundary ( the ratio of specific entropy values above 
and below the direct transition are presented in Fig. 2a). However, it was shown EJ that the 
value {S / Nb)qgp can be corrected by the modification: Bqqp BQGp[^i,T). Indeed, according 
to general thermodynamic relations entropy density s and baryon number density n in QGP are 
defined from: 


nQGp{f^,T) 


9pQGp{^J',T) 

dp 


— T) 


dBQGp{p,T) 

dp 


(5) 


sqgp{p-,T) 


dpQGp{fJ;T,) 

dT 


— ■sqgp(a^> 


dBQGp{p,T) 

dT 


( 6 ) 


where zero superscripts indicate corresponding values calculated for constant B. Other thermody¬ 
namic functions in QGP (e, p, enthalpy w) remain undependent on H’s derivatives. 

Thus to restore conservation of specific entropy it is possible to determine Bqgp{p, T) from 
the differential equation: 


®qgp(A''i 


dBQap{fj.,T) 

dT 

dBQGp{ti.,T) 

diT 


sh{p,T) 

nH{p,T)' 


(7) 


Since S/Nb value is not defined at the points p = 0 and T = 0 it seems natural to fix an integrating 
constant: 


Bqgp{0,T) = Bqgp{p,Q) = Bl^Qp. 

Let us stress that the function Bqgp{p,T) satisfying eq. (7) provides equal behaviour of 
specific entropy functions in both phases everywhere in p—T plane, thus, in particular, conservation 
of the specific entropy at the phase transitioxL-boundary. 

This equation has been solved in Refta; it was shown that the obtained Bqgp{p,T) does 
not change phase diagram considerably. 

However, there still remain several questions concerning this procedure. In particular, 
Bqgp{p,T), as a solution of the eq. (7), is defined in the open region va p — T plane, i.e. for any 
temperature T > Tc{p). Since it depends on the value of specific entropy in H phase, it means 
that EOS of QGP (in particular, the bag constant) should store information on the H-phase EOS 
everywhere including far high-T limit. This looks suspicious since Bqgp represents pressure of the 
physical vacuum and has nothing in common with a particular hadronic system. 

This question becomes serious at relatively large p > 0.9 GeV/fm^ where hadron interactions 
become decisive in the H phaseS. One has to use again phenomenological models for description 
of those interactions (eg., HCM and MFA) which give differing results for the {S/Mb)h value 
(see Fig. 2a). Correspondingly, Bqgp{p,T) modified according to the procedure described above 
should depend on description of nucleon interactions in the H phase. 

Besides, there appears some uncertainty concerning fulfillment of Gibbs relation 

e + p — pn = sT, (8) 

which is to be satisfied in an equilibrium system. Note that this imposes additional constraint on 
B’s derivatives in the case of modified Bqgp(p,T). Indeed, combining (8) with (5),(6) one gets 
the condition: 

.. dB{p,T) _ _rpdB(j^ 

P dfj. ~ dT ^ 

which, seemingly, is not satisfied within the procedure used in 111. 

Another problem concerns dynamical evolution of the system in question. The analysis 
presented refers to the stationary systems. Being related to modern experiments on high energy 
heavy ion collisions, it is to be included in hydrodynamical model of the system evolution (for 
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qualitative analysis the simplest Bjorken scaling solution E3 is used). However the condition (7) 
does not provide equilibrium character of first order phase transition in dynamically evolving 
systems. Such transitions are to proceed through the mixed phase state where not only entropy 
and baryon number are to be conserved, as for all equilibrium processes, but enthalpy of the system 
as well (the latent heat works on the system expansion). 

To provide conservation of these three thermodynamic variables it is necessary and sufficient 
(see Appendix A, eq.(27)) to fulfill the condition: 

WQGpjTc, m /3) _ SQGpjTc, ti/3) _ nQGpjTc, p./3) 
wh(Tc,p) sniTcfJ,) nniTcn) 

Thus, both BQQp{yL,T) derivatives are to be defined according to (9), at least, at the transition 
boundary TSu). Then fixing integrating constant in both ending points, fi — 0 and T = 0, as it 
was done int®l, makes the problem over-defined. 


4 Entropy and baryon number conservation in 
DPTM 

The same problem arises also for transitions considered within DPTM, i.e. for deconfine¬ 
ment and the chiral transition. The ratio S/Nb turns out to be discontinuous when crossing 
both transition boundaries. It is illustrated in Fig.2b, where ’’jumps” of S/Nb ratio at the cor¬ 
responding transition boundaries are presented as calculated, both, in the HCM and the MFA 
phenomenological models describing jmcleon interactions in the H phase. 

Following the same way as inE3, we try to reconsider EOS for deconfined phase. However, 
within DPTM it seems natural and reasonable to modify Q-phase EOS alone making the bag 
pressure parameter Bq /i and T dependent, Bq{pl,T), with QGP bag parameter Bqqp remaining 
constant. 

The reasons are as follows: 

i) . Q-phase is the intermediate phase between H and QGP, thus Q phase EOS can serve as a tool 
for compensation of defects of other phases EOS (remaining the same as earlier) which may occur 
invalid close to the transition boundaries. 

ii) . Bq is chosen even more arbitrarily than Bqq p ■ the last one is to coincide with the QCD vac¬ 

uum pressure estimated usually as 0.5 GeV/fm^, while the Bq is to be closed to the bag pressure 
within the MIT-bag model, and thus varies within the interval: Bq « 50 100 MeV/fm^. 

iii) . The discontinuity of S/Ms is smaller when crossing the deconfinement and the chiral bound¬ 
aries than that for direct transition, (see Eig. 2), thus it is easier to modify Q-phase EOS only. 

iv) . Bq is defined for (and has a physical sense) within closed region of phase space where the Q 
phase can exist {in — T plane: the region bounded by Td{ti) and Tch{p) curves and /r=0 axis). 
Thus there arises no problem with securing proper Bq behaviour at high-T limit, as it appeared 
for BQGp{ti,T). 

v) . modified Bq{pl,T) would also depend on hadron interaction description, and this dependence 
becomes considerable in high-/r limit. However it is quite natural for Bq to be model dependent: 
EOS of the intermediate phase has to vary in accordance with 77-phase EOS variations (in distinc¬ 
tion to model dependence of BQGp{tL,T) in SPTM which seems to be unnatural). 

The modification in question, Bq{pl,T), results in change of entropy and baryon number 
density within Q phase: 


nQ{p.,T) 


dpQ{ti,T) _ 0 ^ dBQ{fi,T) 

dpi dpi 


( 10 ) 




dpQ{ti,T,) _ ^0 dBQ{p,,T) 


( 11 ) 


where zero superscripts indicate corresponding values calculated for constant Bq. Other ther¬ 
modynamic variables in Q-phase (eg, pg, wq) remain unchanged, i.e. do not depend on 77g’s 
derivatives. 
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To provide conservation of specific entropy when crossing transition boundaries and proper 
behavior during the equilibrium mixed phase evolution one needs to fulfill the following conditions 
(see Appendix, eq.(28)): 


WQ{Td,n/3) ^ SQ{Td,n/3) ^ nQ{Td,fi/3) 
WHiTd,n) SH(Td,n) nHiTd,fi) 


WQGpjTch, ^J'/3) _ SQGpjTch, m/3) _ nQGpjTch, ^J■/3) 
WQ{Tch,y^/3) SQ{Tch,y^/3) nQ{Tch,y^/3) 


valid at Q-phase boundaries Td[^i) and Tctifj)- Combining together (10)-(13) we get certain con¬ 
straints on Bq{^,T) function instead of differential equation similar to (7). Namely: 


• at the deconfinement boundary, T(/i) = Td{^j)'. 

i-^)T=TM = ("Q - nH — )T=TMl 

(9Bq. .0 , 

- sh — )t=tm^ 


(14) 

(15) 


at the chiral restoration boundary, T{fi) = TchifJ^) 


y-^)T=TMp) = - riQGP 


.OBq q Wq 

y-^)T=Tch(p) - («Q “ SQGP 


(16) 

(17) 


• at ^ = 0: 

Bq{0,T) = Bq = const (18) 

• inside the Q phase region there is no special constraints on Bq (besides the Gibbs relation 
(8) which is valid for any T and ^), and solution of the problem is not unique, thus the 
function Bq{pl,T) can be chosen rather arbitrarily. 


In classical physics similar problems arise, eg., when simulating soap films stretched on some 
hard contoured (so (piled two-dimensional Plateau problem). Variation methods for such problems 
are well elaborated c3. We have used a numerical procedure providing function Bq{^,T) which 
belongs to the so called minimal surface class, securing local minimum of the surface functional: 


J dfidT^l + a^i?Q(/i,T)2 + 5 ^Bq(^,T)2 


(19) 


with given boundary conditions (14)-(18). Fulfillment of the Gibbs relation (8) has been tested. 

Note that this procedure fails near T — 0 region where the accuracy of calculation becomes 
worse; this case needs special investigation. However we are interested mainly in high-T region 
because this very case can be related to modern experiments on heavy ion ultrarelativistic collisions 
where formation of deconfined phase(s) seems to be rather probable. 


4.1 Results of numerical calculations 

The procedure described has been fulfilled within HGM and MFA models for the following 
parameter values: 

Bqgp = 0.5 GeV/fm3, = 70 MeV/fm^. 
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Fig.3a represents the crossection of the obtained BQ{fi,T) surface for HCM by the decon¬ 
finement transition boundary T = Td{p,) (solid line) and the chiral transition boundary T = TchifJ-) 
(dashed line). The same curves for MFA model are presented in Fig.3b. It is seen that Bq{ii,T) 
remains practically unchanged for small /i values /i < 0.6Gey. For larger p,, Bq{^,T) value 
varies within rather broad interval (relatively to its average value): 50 MeV/fm^ < Bq{^, T) < 90 
MeV/fm^. However the interval of Bq variation is small as compared to difference between Bqcp 
and Bq. Moreover, phase diagrams before and after Bq modification (see Fig.4 ) do not differ 
significantly. The region of Q phase existence remains rather broad with varying Bq as well. 

The modification procedure described above influences mainly the baryon number density 
in the Q phase (see Fig.5), corresponding entropy density corrections satisfy Gibbs relation and 
are relatively small (everywhere except low T region). 

It deserves stressing that BQ{yL,T) differs for HCM and MFA models reflecting intrinsic in¬ 
correctness of the models themselves: Q phase does play its role of intermediate state compensating 
entire defects of H phase EOS. The boundary of chiral transition does not depend practically on 
peculiarities of nucleon interaction in H phase; this means that EOS of QGP does not remember 
the iJ-phase interactions. This seems to be quite reasonable. 


5 Summary and discussion 

It has been shown that within DPTM the modification Bq BQ{yL,T) enables to pro¬ 
vide proper behaviour of thermodynamic functions for reversible equilibrium phase transition, in 
particular, the mixed phase scenario. Phase diagram of three-phases matter is not practically 
changed. The correction concerns mainly baryon number density inside the intermediate Q-phase 
(its changes are not essential). The main result of DPTM - existence of a broad corridor of Q 
phase in fj, — T plane - remains entirely valid. 

Note that within DPTM modifications are necessary and considerable only for sufficiently 
large chemical potential, yi > 0.6-;-0.8 GeV. Eor small /i (most interesting for experimental data 
analysis) the equilibrium character of deconfinement phase transition is almost automatically saved, 
and the change of the chiral transition boundary is negligible. 

Let us stress that EOS of QGP remains unchanged within the method used, thus the problem 
of the direct transition QGP H irreversibility remains as well. But within DPTM the transition 
is to proceed through the intermediate Q phase so that direct transition H QGP should not 
occur normally. However there still remains the possibility for QGP to overcool (too fast) below 
the critical temperature T^h, then the nonequilibrium and thus irreversible phase transition should 
occur. 

It deserves mentioning that the interest to the problem of specific entropy discontinuity has 
been inspired mainly by recent experimental data on heavy ion high energy collision E3 reporting 
a large value of for hadrons resulting from the collision. In this very connection it has 

been pointed out that experimental data does agree with the value typical for QGP and could 
appear in experiment due to abrupt nonequilibrium phase transition QGP —>■ H. However, it 
should be stressed that the ratio S/Nb calculated within DPTM for the resulting hadron gas (see 
Fig.6) is much higher than that for direct transition (and almost the same as in intermediate Q 
phase). It is connected with relatively low temperature of H-phase formation, Td — 140 MeV. Thus 
experimental data could be as well described by equilibrium transition Q —> H instead of abrupt 
direct transition. 

In conclusion-let us point out another problem concerning EOS uncertainties. There were 
put out arguments E3 based on theoretical analysis of effective Lagrangian that hadron (and valon 
as well) masses are to decrease for temperature increasing (since the mass of any hadron is believed 
to be proportional to the quark condensate to the order 1/3), with the fastest decrease (down to 
zero) being close to critical temperature of direct transition in SPTM. Actually this decrease means 
that the deconfined phase with broken chiral symmetry transforms smoothly into QGP without the 
(at least, first order) phase transition. In accordance with such approach DPTM is to be modified 
in such a way that change of valon masses and corresponding change of bag pressure parameters 
(which have to be connected with quark condensate and gluon condensate) are taken into account. 
It seems natural that in this case there occurs the only phase transition, the deconfinement one, 
H ^ Q , while transformation of the Q phase into QGP proceeds smoothly, with constituent 


quark masse approaching zero. Then the problem of the deconfinement transition irreversibility 
remains actual and needs special analysis. 
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Appendix 

Hydrodynamical equations describing the system evolution are to be written for 
energy- momentum tensor; they express the local conservation laws. Neglecting dissipative 
effects, they are: 

= 0, = [e(x) + p{x)]u^ {x)u^ {x) +p{x)g^'^, (20) 

where e is the energy density, p the pressure in the proper system of fluid element, 
the metric tensor (with g^^ = — 1) and is the four velocity of the local flow, x being 
the coordinate of the fluid element in the 4-dimension space. Conservation of the entropy 
and baryon number requires for: 

d^iu'^s) = 0 , ( 21 ) 

dy{v!^n) = 0 , ( 22 ) 

where s and n being respectively the entropy density and baryon number density in the 
proper system of fluid element. In the case of one-dimensional scaling solutionll commonly 
used for qualitative description of evolution of the matter created in heavy ion collisions, 
local conservation of energy-momentum tensor takes the forma (r being the proper time 
of fluid element): 

de _ dr 
e + p T 

Until p (as well as T) is constant during the mixed phase the last equation presents enthalpy 
w = e+p conservation law. Local conservation of entropy and baryon number is expressed 
in a similar form: 

ds dr dn dr 

s T ' n T 

Thus to provide a possibility for an equilibrium transition process through the mixed 
phase state the following conditions are required: 

dw ds dn dr 

w s n T 





In the case of direct transition the mixed phase state is described as a mixture of H 
and QGP fractions taken at constant temperature Tc and given chemical potential p: 

w = wh(Tc, /i)A(r) -b wqgp{Tc, p/3){l - A(r)), (26) 

A(r) being the share of the H-phase admixture which changes during the mixed phase 
expansion from zero (for pure QGP phase) to the unity (for pure hadronic matter). The 
same representation is taken for all additive thermodynamic variables (s, n, e). Thus, (25) 
and (26) combining together lead to the following requirement: 

wqgp(Tc,p/3) ^ SQGpjTcp/S) ^ UQGpjTc, p/3) 
wh{Tc,p) sh{Tc,p) nniTcp) 
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In general case of any first order equilibrium phase transition at some transition 
temperature Ttr corresponding requirement takes the form; 

W+{Ttr,IJ,) ^ S+{Ttr,n) ^ n+(Ttr,fJ,) 

W-{Ttr,n) S-{Ttr,n) n-(Ttr,n)' 

where (+) and (-) refer to thermodynamic variables above and below the transition bound¬ 
ary Ttr in). 

Let us stress that for equilibrium mixed phase scenario it is necessary and sufficient 
to fulfill the conditions (28) only at transition boundaries^ but not for any T and /a. 
Besides, for = 0 (when n=0) and T=0 (when s=0) cases these conditions are satisfied 
automatically due to common thermodynamic Gibbs relation (8). 

It deserves stressing (and can be easily shown) that fulfillment of the conditions (28) 
provides automatically securing of the Gibbs relation (8) in one of the neighboring phases 
if it is valid for the another one. 
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Figure Captions 

Figure 1. Phase diagram in — T plane for nucleon interaction description within HCM (solid) 
and MFA (solid with symbols). Dashed lines correspond to the direct transition in SPTM. 

Figure 2a. Ratio of specific entropy values above (+) and below (-) the direct transition, 
(iS/A/g)’*" : (5 /A/b)“, as a function of baryonic chemical potential /r along the transition boundary. 
Calculated for HCM (solid) and MFA (solid-symbols). 

Figure 2b. The same ratio as in Fig. 2a at the deconfinement (solid) and chiral (dashed) tran¬ 
sition boundaries; curves for MFA model are indicate by symbols. Note the different ordinate scales. 

Figure 3a. Crossection of the obtained Bq{ij,,T) surface by the deconfinement transition bound¬ 
ary T = 7d(/i) (solid line) and the chiral transition boundary T = Tch{ij) (dashed line) as a function 
of baryonic chemical potential fi. Calculated within HCM. 

Figure 3b. The same as in Fig. 3a calculated within MFA. 

Figure 4a. Phase diagram within HCM for Bq=70 MeV/fm^ =Const (dashed) and modified 
BQ{fj,,T) (solid). Short-dashed lines correspond to the limiting values of modified BQ{fj,,T): 
Bq=50 MeV/fm^ and Bq=90 MeV/fm^ for deconfinement and chiral transitions respectively. 
Figure 4b. The same as in Fig. 4a for MFA. 

Figure 5a. Baryon number density in Q phase ng as a function of baryonic chemical potential /i 
at deconfinement (d) and chiral (ch) transition boundaries calculated for HQ=Const (dashed) and 
modified Bq{ij,,T) (solid) for HCM. 

Figure 5b. The same as in Fig. 5a for MFA. 

Figure 6. Entropy per baryon in H phase near the deconfinement in DPTM (solid) and direct in 
SPTM (dashed) transition boundaries as function of quark fugacity A = exp{fXq/T). Short-dashed 
line corresponds to the same ratio in QGP phase. Calculated within HCM. 
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